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Abstract 

In the framework of the large extra dimensions (LED) model, we investigate the effects 
induced by the Kaluza-Klcin gravitons up to the QCD ncxt-to-leading order (NLO) on the 
V7-pair production followed by a subsequcntial W-decay at the CERN LHC. We depict the 
regions in the C — Ms parameter space where the LED effect can and cannot be observed 

from the analyses of the pp —> W + W~ + X and pp — > W + W~ — > W ± l T +X processes. 
We find that the ability of probing the LED effects can be improved by taking the cutoffs 
for the invariant mass of the VU-pair and the transverse momentum of the final lepton. Our 
results demonstrate that the NLO QCD corrections to observables are significant, and do 
not show any improvement for the renormalization/factorization scale uncertainty on the 
QCD NLO corrected cross section, because the leading-order result underestimates the scale 
dependence. 



PACS: ll.10.Kk, 12.38.Bx, 14.70.Fm 
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I. Introduction 



Motivated by the theoretical problems in the standard model (SM), many extended models 
beyond the SM have been established. Among them the large extra dimensions (LED) model 
proposed by Arkani-Hamed, Dimopoulos, and Dvali in Ref.pQ may be one of the promising 
models which can solve the long-standing mass hierarchy problem. This model used the idea 
of extra dimensions to bring gravity effects from the Plank scale down to the electroweak scale. 
In the LED model, the spacetime dimension is D = 4 + 5 with 5 being the dimension of extra 
space, where the gravity and gauge interactions are unified at one fundamental scale M$ ~ TeV 
(the order of the electroweak scale). The graviton propagates in the D-dimensional spacetime, 
while the SM particles exist only in the usual (3 + l)-dimensions. 

Taking into account of the bad behavior of quantum gravity in the ultraviolet (UV) region, 
it is expedient to construct a low-energy effective theory to describe the gravity-gauge-matter 
system in the current (3+l)-dimensional spacetime. In the phenomenological sense, this can 
be achieved through the Kaluza-Klein (KK) reduction in the brane- world scenario [2J. After 
applying this treatment to the LED model, a D-dimensional massless graviton can be perceived 
as a tower of massive KK modes propagating in the (3+l)-dimensional spacetime. It turns out 
that the weakness of gravitational coupling to the SM particles, suppressed by Mp (the reduced 
Planck scale Mp = "^|r=)> can be compensated by summing over numerous KK states. This 
scenario can result in distinct effects at the high-energy colliders [3]. Up to now, many studies 
on the virtual KK graviton effects up to the QCD next-to-leading order (NLO) in the LED 
model have emerged. These include the processes of fermion-pair, multijet, and vector-boson- 
pair production [H El E]- In Ref.[7] the CMS Collaboration has performed a search for LED 
in the diphoton final state events at the \fs = 7 TeV LHC with an integrated luminosity of 
36 pb~ l . They set lower limits on the cutoff scale Ms in the range 1.6 — 2.3 TeV at the 95% 
confidence level. The dijet angular distribution results from the CMS and ATLAS experiments 
appeared in Ref.[8] and provide even stronger limits on Ms, i.e., Ms > 3.4 TeV (CMS) and 
Ms > 3.2 TeV (ATLAS). Recently, the production of a VF-pair at hadronic colliders in the LED 
model has been studied up to the QCD NLO by Neelima Agarwal, et al [Sj. 
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In this paper, we revisit the NLO QCD corrections to the PF-pair production process at the 
LHC in the framework of the LED model, and improve upon the results of Ref.[S] by including 
the effects of top-quark mass and the contribution from the 66-fusion channel. We provide the 
LED effect discovery and exclusion regions, the kinematical distributions up to NLO in QCD 
by taking into account the subsequential Pl^-boson leptonic decay. The rest of the paper is 
organized as follows. In Sec. II, we briefly go into the related Feynman rules in the LED model. 
In Sec. Ill, the leading-order(LO) cross section for the pp — > W + W~ + X process is described. 
In Sec. IV, we calculate the NLO QCD corrections. In Sec. V, we present the numerical results 
for the LO and NLO QCD corrected integrated cross section for the PF-pair production process 
and the distributions of final PF-boson decay products. Finally, a short summary is given. 

II. Related theories 

The LED model consists of the pure gravity sector and the SM sector. In this model the manifold, 
in which gravity propagates, is not the ordinary four-dimensional spacetime manifold M 4 , but 
R 4 x M, where Ai is a compact manifold of dimension 5. For simplicity, one can tentatively 
assume that M is a 5-torus with radius R and volume Vs = (2ttR) s without loss of physical 
significance. 

In our work we use the de Donder gauge. The Feynman rules for the propagator of the spin- 
2 KK graviton and the relevant vertices which we use are listed below. There Gg K , V, W^, 
A aiM , and rf represent the fields of the graviton, quark, PF-boson, gluon, and SU(3) ghost, 
respectively. 

• spin-2 KK graviton propagator after summation over KK states : 



D - 2 
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(2.1) 



• ^kk^) - iPifa) ~ ip(ki) vertex : 

- i-=L- [r(h + k 2 )» + y (fci + k 2 r - 2rT(h + h~ 2m^)] 



(2.2) 



• G^ K (k 4 ) - ij(h) - i>(k 2 ) - A a P(k 3 ) vertex : 




27V") T 



(2.3) 
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G^ K (k 3 ) - A a P(k±) - A ba {k 2 ) vertex : 



M P 



"3 



G^ K (k 3 ) - W + P{h) - W~ a (k 2 ) vertex : 



M P 



BP upa m 2 w + (C^ CTT/3 - C^ aT )k lT k 2 p + jE^fa, k 2 



G^ K (k 4 ) - A a P{h) - A b(7 (k 2 ) - A cX {k 3 ) vertex : 



W P 9sf 



abc 



(fei - k 3 ) T CP UTapx + (k 2 - k x ) T CP vapTX + {k 3 - k 2 ) T CP uXaTp 



Gg K (k 5 ) - A a P(h) - A b °{k 2 ) - A cX {k 3 ) - A dS {k A ) vetex 



^ _ g2^jeac jebd j^pvpaXS _|_ yeafe j:ecdj-^p,upXaS _|_ jead jebc jjpi/pa 

M p 



5X] 



G KK( k s) ~ V a (h) - r, b (k 2 ) vertex : 



i 5 ^B a ^k la k w 
M p 



G kk(^) - V a (h) - r, b (k 2 ) - A c P(k 3 ) vertex : 



(2.4) 



(2.5) 



(2.6) 



(2.7) 



(2.8) 



M P 



gjabczappv^ 



(2.9) 



where g s is the strong coupling constant, T a and f abc are SU(3) generators and structure con- 
stants, D = n + 5, n = 4 — 2e, M p is the reduced Planck mass, 03 and £ are SU(3) and charged 
SU(2) gauge fixing parameters, and D(s) can be expressed as [2] 



D(s) 



T(5/2) (4vr) 5 / 2 



7T + 2i/(A/ v / i) 



and 



r-A/v^ 5-1 

I(A/V~s) = P I dy 



(2.10) 



(2.11) 



/o 1 - 

The integral I(A/^/s) contains an ultraviolet cutoff A on the KK modes [21 E]- In this work we 
set it to be the fundamental scale Ms- It should be understood that a point y = 1 has been 
removed from the integration path. Besides, all the momenta are assumed to be incoming to 



1 

(1) (2) (3) 

Figure 1: The tree-level Feynman diagrams for the partonic processes qq — > W + W~ in the LED 
model. (1) and (2) are the SM-like diagrams, where q represents the u-, d-, o, s- and 6-quark. 
(3) is the extra diagram with KK graviton exchange. 



the vertices, except that the fermionic momenta are set to be along the fermion flow directions. 
The coefficients A pu , B pvOL P , C pap f ia ' 3 ', D pu P aXS , and E^P a (k 1 ,k 2 ) are expressed as 

D ^paXS = ^u^pa^XS _ rf&rf^ + (jfP^^ + rf x rf u r] p& - rf p rf a r] X& - rf x rf & r) pa + (/x O u)), 
E pupa (k u k 2 ) = Tf v (k{kl + k p 2 k% + k p kl) - [if a k p k{ + jf p k^kl + (/x <-> v)\ . 

We code programmatically the related Feynman rules in the FeynArts 3.5 package |10j to 
generate the Feynman diagrams and the relevant amplitudes. The FormCalc 5.4 [11] package is 
implemented subsequently to simplify the amplitudes. 

III. LO cross section for pp — > W + W + X 

We treat the up-, down-, charm-, strange-, and bottom-quark as massless particles, and adopt the 
five-flavor scheme in the leading order and QCD next-to-leading order calculations. The LO con- 
tribution to the parent process pp — > W + W~ +X includes the quark-antiquark (q = u, d, s, c, b) 
annihilations and the gluon-gluon fusion partonic processes: q(pi) +q{p2) W + {pz) + W~(p4) 
and g{p\) + g(j>2) — > W + {ps) + W~{p^). There pi (i = 1,2,3,4) represent the four-momenta 
of the incoming and outgoing particles, respectively. The corresponding Feynman diagrams are 
shown in Figsfl]andG±l FiguresQJl) and H^2) are the LO SM-like diagrams for partonic process 
qq ->■ W + W~. In FigOJl) the internal wavy line means exchanging 7 or a Z°-b oson. There we 
ignore the diagrams with exchanging Higgs boson, since the initial quarks are all massless. 
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Figure 2: The tree-level Feynman diagram for the partonic process gg — > W + W in the LED 
model. 



We express the tree-level amplitudes for the partonic processes qq — > W + W and gg — > 
W + W~ as 

M% = M% SM + M% LED , M% = M% LED , (3.1) 

where A4®'^ M (q = u, d, c, s, b) is the amplitude contributed by the treedevel SM-like diagrams, 
while M.% ED and M?gg ED are the treedevel amplitudes with KK graviton exchange. Our 
calculations show that the analytical expression of the SM matrix element squares summed 



(averaged) over the final (initial) state spins and colors at the LO, | 2 , for the partonic 

process without massive internal or external quark (i.e., q = u,d,c,s), is the same as that 
presented in Ref.[9]. But for the partonic process bb — > W + W~, there is a t-channel diagram 



with a massive top-quark exchange. Its explicit expression of l-Ad^g 5 *^ 2 is presented below by 
adopting the notations in Ref . [9] . 

lKf M ! 2 = ^ ( A l B lm + AB b 2m + > (3-2) 

where N is the number of colors. The explicit expressions for A\, A^, and can be obtained 



from the Eqs.(8) in Ref.[9j by taking the replacement of u — > b. The kinematic invariants B\ m , 
B\ m , and B^ m can be expressed as 



2 

B\ m = . U 2 . 2 Bf(t,u,s), B b 2m = Bl(t,u,s), B\ m = - B$(t,u,s), (3.3) 
(u — rcitY [u — mi) 



where Bf(t,u,s), B$(t,u,s) and B$(t, u,s) are presented in Eqs.(9)-(12) of Ref.[9]. Then the 
LO cross sections for the unpolarized VF-pair production processes at the partonic level can be 
expressed as 

vjf = -^-^ j dT 2 \M%\\ (ij = uu, dd, cc, ss, bb, gg) , (3.4) 
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where p is the momentum of one initial parton in center-of-mass system (c.m.s.) and is the 
two-body phase space element expressed as 

dr 2 = (27r)V 4 )( Pl + P2 - P3 - PA ) J] jS^r.- ^ 

i=3,4 ^ ' 1 

By convoluting with the parton distribution functions (PDFs) of the colliding protons, 
the LO cross section for the parent process, pp — > W + W~ + X , can be written as 

cc.bb.gg ^ „ 

G i/A (x A ,(if)G j/B (x B ,Vf)aif(\f§) + (*<-> j) , 

ij=uu,dd,ss J 

(3.6) 

where G^/p {i = g,q,q) represent the PDFs of parton i in proton P, /if is the factorization 
scale, \/I = xaxbVs, %A and xb describe the momentum fractions of parton (gluon or quark) 
in protons A and B, respectively. 

IV. NLO QCD corrections 

The complete NLO QCD correction to the parent process pp — > W + W~ + X consists of follow- 
ing components. (1) The virtual contribution from the QCD one-loop and the corresponding 
counterterm diagrams to the partonic channels qq -> and gg — > W + W~. (2) The 

contribution of the real gluon emission partonic processes. (3) The contribution of the real 
light- (anti) quark emission partonic processes. And (4) the corresponding contribution of the 
PDF counterterms. There inevitably exist the ultraviolet (UV) and infrared (IR) divergences 
in the NLO calculations, and we adopt the dimensional regularization scheme in n = 4 — 2e 
dimensions to isolate and manipulate these divergences. 
A. Virtual corrections 

The Feynman diagrams for the virtual corrections to the qq — > W + W~ and gg — > W + W~ partonic 
processes are shown in FigjH] and FigJU respectively. In Figs|3[3) and IU(4) the diagrams in- 
volving Yukawa coupling between Higgs boson and top quarks are included, but the diagrams 
involving Yukawa coupling between Higgs boson and massless quarks are excluded due to their 
vanishing contribution. There exist UV and soft/collinear IR singularities in the calculations 
of these one-loop diagrams. To remove the UV divergences, we need only the wave function 
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E 



1 



dxAdxs 



renormalization constants for the quark and gluon fields. We introduce the renormalization 
constants 5Z^ L R for massless quark (q=u,d,c,s,b) fields and 5Za for the gluon field defined as 

i/, q>L>R = (1 + SZ^) 1 /^^ Af = (1 + SZa) 1 /^. (4.1) 



In the modified minimal subtraction (MS) renormalization scheme the renormalization constants 
for the massless quarks are expressed as 

sz *,.l = ~C F (A UV - A IR ), SZ^ R =-^C F (A„ V -A IR ), (4.2) 

sz * = S - r v ' VT ") Avv + £ (! c * - t'" Tf ) Ai «- (4 - 3) 

To remove the UV and IR divergences in the 66-fusion subprocess, we need introduce the coun- 
terterms for the top-quark field and its mass, i.e., 

1>t,L,R = 0- + 8Zf ttLtR ) 1 / 2 il>t, LtR , m° t =m t + 5m t . (4.4) 

We use the on-mass-shell scheme to renormalize the top-quark field and mass. They are expressed 
as 

2 

^ Z 4h,l = -gc F (A c/v + 2A /R + 31n^+4), (4.5) 
5Z^ R = -^CHA c/v + 2A JR + 31n^| + 4), (4-6) 



5m t 3a s . \i 



C F (A uv + \n^ + -), (4.7) 
nit 47r mi o 

In the above equations \i r is the renormalization scale, Cp = |, Ca = 3, Tp = g, n^ v = 6 
corresponds to the six flavor quarks (u, d, c, s, t, b), whereas n^ R = 5 is the number of 
the massless quarks (u, d, s, c, b). Moreover, Auv = ^^T(l + euv){^Y uv and Air = 
^j^r(l + em)(4:ir) eiR refer to the UV and IR divergences, respectively. 

Then the results for the differential cross sections for the qq annihilation and gg fusion par- 
tonic channels are UV finite but soft/collinear IR divergent. The soft /collinear IR singularities 
can be canceled by adding the contributions of the real emission partonic processes and the 
corresponding PDF counterterms. 

B. Real gluon emission 



S 



1 11 '/ •/ 1 




(5) (6) (7) (8) 




(9) (10) 



Figure 3: The QCD one-loop Feynman diagrams for the partonic process qq — > W + W . (l)-(6) 
are the SM-like diagrams. (7)-(10) are the diagrams with KK graviton exchange. 

The real gluon emission contributions are from g(p\) + g{p2) — > W + (ps) + W~{pi) + g{p§) 
and q(pi) + q(p2) — > W + (ps) + W~{p^j + g(p§) partonic processes. The corresponding Feynman 
diagrams are shown in Figj5] and Figj6l respectively. We employ the two cutoff phase space 
slicing (TCPSS) method [12j to calculate the contributions from the real gluon emission partonic 
processes. An arbitrary soft cutoff 5 S is introduced to separate the gluon emission subprocess 
phase space into two regions, soft gluon and hard gluon regions. Furthermore, another cutoff 5 C 
is introduced to decompose the real hard gluon emission phase space region into hard collinear 
(HC) and hard noncollinear (HC) regions. The partonic differential cross section for the real 
gluon emission subprocess can be expressed as 

dd g = da s g +da^ = da s g + daf C + daf c '. (4.8) 
C. Real light-(anti)quark emission 

In addition to the real gluon emission discussed above, there are contributions from the 
massless light-(anti)quark (u,d,c,s,b,u,d,c,s,b) emission partonic processes. In the five-flavor 
scheme the massless light-quark q involves u-, d-, o, s-, 6-quarks. Considering the fact that 
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(1) 



(2) 



(3) 



(4) 




(9) 



(10) 



(ID 



(12) 




(13) 



(14) 



(15) 



(16) 




(17) 



(18) 



(19) 



(20) 



Figure 4: The QCD one-loop Feynman diagrams for partonic process gg — > W + W~ . (l)-(7) are 
the SM-like diagrams. (8)-(20) are the diagrams with KK graviton exchange. In all diagrams q m 
represents u-, d-, c-, s-, b- and t-quark except the diagrams in FigsH](3) and Fig2j4) involving 
the coupling between Higgs boson and top quarks, where q m denotes only top quark. 
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the final (anti)bottom-quark can be tagged in experiments and the collinear IR singularities 
of the real (anti)bottom-quark emission subprocesses are completely canceled by those of the 
corresponding PDF counter terms, we do not include the contributions of the bottom and 
antibottom emissions, and adopt the five-flavor PDFs [13] . We depict the Feynman diagrams 
for the partonic processes qg —> W + W~ + q and qg —> W + W~ + q in FigJT] These partonic 
processes contain only the initial state collinear singularities. By using the TCPSS method 
described above, we can split the phase space into collinear and noncollinear regions. The 
differential cross sections for the partonic processes qg —> W + W~q and qg — > W + W~q can then 
be expressed as 

da q{q) (q(q)g^W + W-q{q)) = da c q[q) + da^ (4.9) 

where da^ and d&^ are finite and can be evaluated by using the general Monte Carlo method. 
D. NLO QCD corrections to the pp — > W + W~ + X process 

Combining the renormalized virtual corrections and the real gluon/light-(anti)quark emission 
contributions, the partonic cross sections still contain the collinear divergence, which can be 
absorbed into the redefinition of the PDFs at the NLO according to the mass factorization [14] . 
We find that after the summation of all the NLO QCD corrections, the soft and collinear IR 
divergences vanish. We can see from above discussion that the final total 0(a s ) corrections 
consist of the two-body term and the three-body term a^ 3 \ The total cross section up to 
the QCD NLO is expressed as 

(J NLO = ctlo + Actjvlo = <?LO + cr {2) + <t (3) . (4.10) 

It is UV finite, IR safe, and cutoff 5 C /5 S independent |12[ 115], which will further be checked in 
the numerical evaluations. 

V. Numerical results and discussions 

In this section, we present the numerical results of the integrated cross sections and the kinematic 
distributions of the final particles for the pp — > W + W~ + X process in both the SM and LED 
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(1) 



(2) 



(3) 




(7) (8) 



Figure 7: The tree- level Feynman diagrams for the real light- (anti) quark emission subprocesses 
q(q)g — > W + W~ q(q) . (l)-(5) are the SM-like diagrams. (6)-(9) are the diagrams with KK 
graviton exchange. 
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model up to the QCD NLO. In order to verify the correctness of our numerical calculations, we 
made the following checks: 

(i) In Table HJ we list our numerical results of the LO and NLO QCD corrected integrated 
cross sections in the SM for the pp — > W + W~ +X process by taking the input parameters, 
PDFs, and event selection criterion from Table 4 of Ref.[16]. It shows that our LO and 
NLO QCD corrected cross sections in the SM are in good agreement with those in Ref.[16| 
within the integration errors. 



LHC 


Ref.[16j 


FeynArts 


CompHEP 


Born[pb] 


86.7 


86.711(6) 


86.7(1) 


NLO QCD[pb] 


127.8 


127.7(1) 





Table 1: The LO and NLO QCD corrected cross sections for the pp —> W + W + X process in 
the SM by taking the same input parameters and event selection criterion as those in Ref.[16j. 



(ii) The UV and IR safeties are verified numerically after combining all the contributions at 
the QCD NLO. 

(in) We calculate the NLO QCD corrections to integrated cross section for the pp —> uu — ?■ 
W + W~ + X process as functions of the cutoff 5 S at the ■s/s = 14 TeV LHC in the LED 
model, where we take fj,f = \x r = fio = mw, Ms = 3.5 TeV, 5 = 3 and 5 C = 5 s /50. Some of 
the results are listed in Table [2j It is shown clearly that the NLO QCD correction (ActJ^q) 
does not depend on the arbitrarily chosen values of 5 S and 5 C within the calculation errors. 
In the further numerical calculations, we fix 5 S = 1(T 3 and 5 C = 6 s /50. 

(iv) We calculate the SM LO VF-pair invariant mass distribution {daf^ / dMy/w) for the pp — > 
W + W~ + X process with the same input parameters, PDFs and event selection criterion 
as those used in Ref.0. The numerical results, which are obtained by using FeynArts and 
CompHEP packages separately, are coincident with each other. 

(v) For further comparison with the previous work of N. Agarwal, et al, we recalculate the LO 
and NLO QCD corrected distributions of W^-pair invariant mass in both the SM and LED 
model (da^' Lt,L> /dMww, da^LO /dM ww ) for the pp ->■ W + W~ + X process at the 
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?5 

2 x lO -3 


13.19(3) 


?5 — 

1 x 10~ 3 


13.20(3) 


7 x 10" 4 


13.21(3) 


4 x 10~ 4 


13.22(5) 


2 x 10~ 4 


13.24(5) 


1 x 10" 4 


13.26(6) 


7 x 1(T 5 


13.25(6) 


4 x 10~ b 


13.27(6) 


2 x 1(T 5 


13.26(7) 



Table 2: The dependence of the NLO QCD correction to the integrated cross section for the 
process pp — > uu — > W + W~ + X at the ^fs = 14 TeV LHC in the LED model, where we set 
fj,j = fi r = fio = m\Y, Ms = 3.5 TeV, 5 = 3 and 5 C = <5 s /50. 

y/s = 14 TeV LHC, where we set all the quarks being massless except mj = 172.0 GeV, 
and take the PDFs and event selection criterion from Ref.[9]. We plot our LO and NLO 
QCD corrected results in the SM in FiglS^a), and the results in the LED model in Figj8](b). 
In these two figures we depict also the corresponding curves from N. Agarwal's paper [9] 
for comparison. We can see that there exist obvious discrepancies between ours and the 
corresponding N. Agarwal results, especially in the large Mww region. One of the reasons 
for the disagreement is because we have included the effects of top-quark mass in our 
calculations. From Figsj8|a,b) we can see the K-factors of the QCD corrections increase 
to large numbers of K = 5.29 and K = 2.33 separately, when VF-pair invariant mass 
approaches My/w = 1300 GeV . This occurs because the K-factors in FigslU(a,b) are the 
results with only constraint on W-bosons (ly^l < 2.5 ) [9j. In this case the differential 
cross section, daf/^Q/dMww, receives a large contribution from the hard jet emission 
corrections (W" 1 "!^ - + jet). For example, in Figj8|a) at M\yw = 1300 GeV the K-factor 
contributed by hard jet emission processes can reach 3.59. 

In the following calculations we take one-loop and two-loop running a s in the LO and QCD 
NLO calculations, respectively [17] • The QCD parameters are taken as = 165 MeV, 

\MS _ 226 MeV, the number of the active flavors is rif = 5, the Cabibbo-Kobayashi-Maskawa 
matrix is set as a unit matrix, and the colliding c.m.s. energy is taken as y/s = 7 TeV and 
y/s = 14 TeV for the early and future LHC. To satisfy the unitary constraint, we adopt the 
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Figure 8: Invariant mass (M\yw) distributions at the LO and NLO for the pp — > W + W~ + X 
process at the ^fs = 14 TeV LHC. There the input parameters, PDFs, and event selection 
criterion are taken from Ref.0 (where Ms = 2 TeV, 5 = 3, and fj, r = fj,f = My/w)- For 
comparison with previous work, we depict also the corresponding curves from N. Agarwal's 
paper [9] in both panels, (a) The distributions in the SM. (b) The distributions in the LED 
model. 

cut t/S < Ms for the whole phase space. We assume mu = 120 GeV and the renormalization 
and factorization scales to be equal (/i r = fif = fi), and we define fiQ = m\\r. We use the 
CTEQ6L1 and CTEQ6M PDFs [MIS] in the LO and QCD NLO calculations, respectively. The 
other related input parameters are taken from [17j : a~ 1 (mz) = 127.916, mw = 80.399 GeV, 
mz = 91.1876 GeV and mt = 172.0 GeV . As we know that the constraints on the final 
particles are necessary in realistic experimental event collections, and the theoretical calculation 
should keep the perturbative convergence. We adopt the following event selection constraints 
additionally. (1) For the real emission contributions, we accept the events which satisfy the 
condition that the jet pseudorapidity \yj e t\ > 2.5 or the jet transverse momentum p^f* < 50 GeV. 
(2) The VK-pair invariant mass is restricted in the range of Mw\y > 400 GeV. 

In FigsHIa,b), the upper panels show the scale (/u) dependence of the LO and the NLO 
QCD corrected cross sections in the SM and LED model at the y/s = 7 GeV and y/s = 14 TeV 
LHC separately, and the corresponding if-factors (K(fj,) = ^^rjjr ) are illustrated in the lower 
panels. There we take Ms = 3.5 TeV and 5 = 3. The scale-dependent .fT-factor in the LED 
model varies from 1.18 (1.53) to 1.19 (1.11) when [i goes from 0.5//o to at the ^fs = 7 TeV 
LHC (the ^/s = 14 TeV LHC). We see from these upper panels that the NLO QCD corrections 
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in the SM and LED model do not reduce the factorization/renormalization scale uncertainty, 
especially at the y/s = 14 TeV LHC. That is because the LO result underestimates the scale 
dependence due to the LO contribution being from pure electroweak partonic processes. And 
we find that the A-factors plotted in the lower figures keep the convergence of the perturbative 
series in the plotted \x range at both the y/s = 7 TeV and y/s = 14 TeV LHC. In further 
calculations we fix /x = mw- 




Figure 9: The scale dependence of the LO and NLO QCD corrected cross sections for the process 
pp — > W + W~ +X in the SM and LED model, and the corresponding if -factor [A(/x) = ^"^ff ] 
with M s = 3.5 TeV and <5 = 3. (a) At the y/s = 7 TeV LHC. (b) At the y/s = 14 TeV LHC. 



In Figs llO( a.b). we depict the LO and NLO QCD corrected cross sections and the corre- 
sponding A-factors for the process pp — > W + W~ +X in the LED model as the functions of the 
fundamental scale Ms, with ^ = mw and the extra space dimension number 5 being 3, 4, and 
5, respectively. From the figures one can find that the more distinct LED effect exhibits with 
the smaller values of Ms and 5. 

The LO and NLO QCD corrected distributions of the VF-pair invariant mass and the cor- 
responding A'-factors (K(M WW ) = dM^w ^ dMww ^ for the P rocess PP ~^ W + W~ + X in the 
SM and LED model at the early and future LHC, are shown in Figslll(a) and (b), separately. 
There the results are for Ms = 3.5 TeV, /i = mw, at a fixed value 3 for the number of extra 
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Ms [TeV] Ms [TeV] 



Figure 10: The LO and NLO QCD corrected cross sections and the corresponding K-factors 
for the process pp — > W + W~ + X in the LED model as functions of Ms with \i = my/ and 
5 = 3,4, 5. (a) At the y/s = 7 TeV LHC. (b) At the y/s = 14 TeV LHC. 

dimensions and obtained by taking the input parameters and the event selection constraints 
mentioned above. As we expected, the LO and NLO QCD corrected differential cross sections 
of the VF-pair invariant mass become less with the increment of Mww- 

It is clear that if the deviation of the cross section from the SM prediction is large enough, 
the LED effect including the NLO QCD corrections can be found. We assume that the LED 
effect can and cannot be observed, only if 



= \°NLO ~ a NLO\ > £ = 5fT I 5 - 1 ] 



and 



a \„LED „SM I ^ 3y£cr NLO _ 

A-CTNLO = WnLO ~ a NLO\ < £ = dcr : K b - 2 ) 

respectively. In FigsJ12fa.b). we present the discovery and exclusion regions in the luminosity- 
fundamental scale space (C — Ms) for the pp — > W + W~ + X process with 5 = 3. FigurefT^Ta) 
is for the = 7 TeV LHC and FigH^b) for the ^ = 14 TeV LHC, where the LED effect 
can and cannot be observed in the dark and gray- region, separately. We list some typical data 
which are read out from FigsJ12fa.b) in Table There the discovery and exclusion fundamental 
scale Ms values at the early and future LHC are presented. It shows that by using the W-boson 
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Figure 11: The LO and NLO QCD corrected distributions of the W-pair invariant mass 
(d(TLo/dMww ; doNLo/dMy/w) and the corresponding if-factors for pp — > W + W~ + X with 
M s = 3.5 TeV, /i = m w and 5 = 3 in the SM and LED model, (a) At the y/s = 7 TeV LHC. 
(b) At the y/s = 14 TeV LHC. 

pair production events we could set an exclusion limit on the cutoff scale Ms to be 1.80 TeV at 
the 95% confidence level at the y/s = 7 TeV LHC with an integrated luminosity of 36 
This is in the Ms lower limit range of 1.6 ~ 2.3 TeV obtained experimentally by the CMS using 
the diphoton final state data samples [7j. 



Luminosity(£) 

Oc) 


y/s = 7 TeV 


x/s = U TeV 


M s [TeV](3a) 


M s [TeV](5a) 


Ms[TeV]{3a) 


M s [TeV](5a) 


100 fb~ L 


3.83 


3.50 


5.69 


5.46 


200 fb~ L 


3.98 


3.74 


5.79 


5.62 


300 fb~ L 


4.17 


3.85 


5.83 


5.70 


36 pfe -1 


1.80 


1.68 


2.03 


1.89 



Table 3: The discovery (Acrtot > 5cr) and exclusion (Aat t < 3<r) LED model fundamental 
scale (Ms) values for the pp — > W + W~ + X process at the early (y/s = 7 TeV) and future 
(y/s = 14 TeV) LHC. 



Now we consider the subsequential leptonic (electron, muon) decay of one of the two W- 
bosons. In collecting the pp — > W + W~ — > W^P t'j +X (I = e, fi) events we do not distinguish 
the leptonic charge. We fix the branching fraction for VF-boson decay (W^ — > 1^ v , i = e,fi) 
as 21.32% j!7j . C = 300 fb" 1 , and take the number of the extra dimensions (5 = 3, the constraints 
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Figure 12: The LED effect discovery area (dark) and the exclusion area (gray) in the C — Ms 
space for the pp ->■ W + W~ +X process, (a) At the y/s = 7 TeV LHC. (b) At the y/s = 14 TeV 
LHC. 

of M\vw > 400 GeV, p l T > p™\ = 100 GeV, and the jet event selection criterion as declared 
above. We show the discovery and exclusion regions in the C — Mg space for the processes 
pp -> W + W~ -> W^fF W +X (I = e,n) in FigsfH^a) and (b) for the y/s = 7 TeV and 
y/s = 14 TeV LHC, respectively, The dark and gray- regions represent the parameter space 
where the LED effect can and cannot be observed separately. Some representative data for the 
discovery and exclusive fundamental scale Ms values at the early (y/s = 7 TeV) and future 
(y/s = 14 TeV) LHC read out from Figs ll3( a.b) are presented in Table HI We can see from the 
table that by using pp — > W + W~ — > W^P v\ +X (I = e, fx) processes with the constraints of 
M ww > 400 GeV, p l T > pf\ = 100 GeV, and our chosen jet event selection criterion, we could 
get exclusion lower limit on Ms as 2.19 TeV at the 95% confidence level at the y/s = 7 TeV 
LHC with an integrated luminosity of 36 (pb)' 1 , which is larger than that obtained by analyzing 
the VF-pair production events as described above. 

We depict the LED discovery and exclusion regions in the Ms — Pjtf space for the processes 
pp -»• W + W~ -»• W ± P y +X (t = e,n) in FigsIH^a) and (b) with 5 = 3, C = 300 fb~ l , 
Mww > 400 GeV, and the branching fraction for VF-boson decays (W T — > £ T v , I = e, fj,) 
as 21.32%, where Fig{TJJa) and FigHHb) are for the y/s~ = 7 TeV and y/s~ = 14 TeV LHC 
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Figure 13: The LED effect discovery area (dark) and exclusion area (gray) in the C — Mg 

space for the pp — > W + W~ — > W^l T V[ +X (£ = e, /x) processes with the constraints of 
M\yw > 400 GeV, p l T > pjtf = 100 GeV, and the jet event selection criterion declared above, 
(a) at the y/s = 7 TeV LHC. (b) at the y/s = 14 TeV LHC. 



Luminosity(zC) 

(£) 


= 7TeV 


y/s = UTeV 


M s [TeV] (3(j) 


M s [TeV](5a) 


M s [TeV](3a) 


M s [TeV](5a) 


100 fb- 1 


4.42 


3.95 


6.69 


6.41 


200 fb~ l 


4.65 


4.29 


6.82 


6.62 


300 fb- 1 


4.74 


4.45 


6.87 


6.71 


36 pb~ L 


2.19 


1.96 


2.98 


2.80 



Table 4: The discovery (Aatot > 5<t) and exclusion (Aa to t < 3<r) fundamental scale (-Ms) 

values for the pp -> -> V/ +X (4 = e,/i) processes in the C — Ms space for the 

pp — > W + W~ — > W ± l T ui +X (£ = e,/i) processes with the constraints of Mww > 400 GeV 
and p l T > pfj = 100 GeV: at the y/s = 7 TeV LHC and at the y/s = 14 TeV LHC. 
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respectively. The dark and gray- regions represent the parameter regions where the LED effect 
can and cannot be observed, separately, with the constraints of prj > Pj^j and the VF-pair 
invariant mass M\y W > 400 GeV. Some representative data are listed in TableOfor the discovery 
and exclusion fundamental scale Ms values with different pjff values at the y/s = 7 TeV and 
y/s = 14 TeV LHC as shown in Figsll4(a.b). We can see that in the case where we fix the 
integral luminosity (e.g. C = 300 we could improve slightly the low limit on Ms if we 

adopt a larger lower cut on lepton transverse momentum (pj^j)- 




Figure 14: The LED effect discovery area (dark) and exclusion area (gray) in the Ms — Pj^f 

space for the pp — >■ W + W~ — >■ W ± l T ^v} +X (£ = e, fj) processes with 5 = 3 and C = 300 fb~ l . 
(a) At the y/s = 7 TeV LHC. (b) At the y/s = 14 TeV LHC. 



p l T cut value 
(Ptj) 


7 TeV 


14 TeV 


M s [TeV](3a) 


M s [TeV](5a) 


M s [TeV](3a) 


M s [TeV](5a) 


50 GeV 


4.61 


4.24 


6.67 


6.60 


100 GeV 


4.74 


4.45 


6.87 


6.71 


150 GeV 


4.92 


4.61 


6.92 


6.80 



Table 5: The discovery (Aatot > 5cr) and exclusion (Aatot < 3cr) LED model fundamental scale 
(Ms) values in the Ms — p™f space for the pp -»• W + W~ -»• W ± P u{ +X (I = e, /z) processes 
with the constraints of M ww > 400 GeV and p l T > pfj. (a) at the y/s = 7 TeV LHC. (b) at 
the y/s = 14 TeV LHC. 
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VI. Summary 

We calculate the NLO QCD corrections to the pp — > W + W~ — > W^P +X process in the 
SM and LED model at the LHC. We investigate the integrated cross sections, the distributions of 
some kinematic variables and how they are affected by radiative corrections. The calculations are 
compared with previous works, and finally the reliable numerical results are obtained. We find 
that the NLO QCD corrected results do not show remarkable reduction of the scale uncertainties 
of the LO cross sections in both the SM and LED model, because the uncertainty of the LO 
cross section is underestimated. The scale-dependent iT-factor is found to be the value from 
1.18 (1.53) to 1.19 (1.11) when /i goes from 0.5^ to 2/i at the sfs = 7 TeV (y/s = 14 TeV) 
LHC, with the constraints of Myyw > 400 GeV and our jet event selection criterion. The 5a 
discovery and 3d exclusion ranges for the LED parameters Ms are also obtained in the NLO 
QCD. The inclusion of the effects of the virtual KK graviton turns out to enhance the differential 
distributions of kinematical observables generally. We conclude that the NLO QCD correction 
to the VF-pair production make it possible to precisely test the TeV quantum gravity in the 
LED scenario at the LHC. 
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